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The plasmon damping rate in scalar field theory is com- 
puted close to the critical temperature. It is shown that the 
divergent result obtained in perturbation theory is a conse- 
quence of neglecting the thermal renormalization of the cou- 
pling. Taking this effect into account, a vanishing damping 
rate is obtained, leading to the critical slowing down of the 
equilibration process. 
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The behaviour of field theory at finite temperature has 
been the subject of intense study in the recent past, mainly 
in connection with the investigation of phase transitions in 
the early universe and of the properties of the quark-gluon 
plasma. Most of the effort has been devoted to the calcula- 
tion of time-independent quantities, like the QCD pressure or 
the electroweak effective potential. Thus, finite temperature 
field theory in thermal equilibrium has been mainly employed 
in this kind of studies. 

More recently, time-dependent processes have been taken in 
consideration, like the generation of the cosmological baryon 
asymmetry or the dynamics of scalar fields in inflationary 
models. 

An example of a dynamical quantity which has received 
much attention in the recent literature is the plasmon damp- 
ing rate, defined as 
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where m(T) is the thermal mass and II/ is the imaginary part 
of the two-point retarded green function. 

7(T) gives the inverse relaxation time of space independent 
fluctuations at the temperature T [jjj, and in A$ 4 /4! theor 
it emerges at two-loops in perturbation theory. In refs 
it was computed as 
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where l qu = \/m(T) — ^35-^ is the Compton wave- 
length, and the temperature T is much larger than any mass 
scale of the T = theory. 

In refs. Q it was shown that the above two-loop result can 
be reproduced in the classical theory provided that the Comp- 
ton wavelength l qu is identified with the classical correlation 
length l c . The above result can be understood realizing that 
7 probes the theory at scales uj — m(T) <C T (if A is pertur- 
batively small) , where the Bose-Einstein distribution function 
is approximated by its classical limit, 
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with p — 1/T. Considering for instance the 1-1 component of 
the propagator in the real-time formalism, 

An = " 2 ^ 5 ( fc2 - ™) (5 + * (M)) , (4) 

we see that when the loop momenta are ko <C T the 'statis- 
tical' part of the imaginary part of the propagator dominates 
over the 'quantum' part, i.e. N 1/2, and the leading or- 
der result can be obtained neglecting the T = quantum 
contributions to the loop corrections. 

The above argument has been employed to motivate the use 
of classical equations of motion in the study of the evolution 
of long wavelength modes in scalar and gauge theories. More 
recently, this approach has been improved by many authors 
including the effect of Hard Thermal Loops in the equations 
of motion for the 'soft' modes (see for instance ||[|). The sep- 
aration between hard and soft modes has been made explicit 
in ref. || by introducing a cutoff A such that m(T) A < T. 

In this letter we want to investigate further the relation 
between thermal and classical field theory, extending the 
computation of the damping rate to temperatures close to 
the critical one. The 1-loop thermal mass is now given by 
m 2 (T) _ + xt 2 /24. In the limit 
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the thermal mass vanishes and the behaviour of the infrared 
(IR) modes is exactly statistical, see eq. (||). 

The divergence of the correlation length at the critical point 
implies that the two- loop/classical expression in eq. (fe|) di- 
verges as well. Physically, this would mean that the lifetime 
of long wavelength fluctuations is getting shorter and shorter 
as the critical temperature is approached, a behaviour oppo- 
site to the critical slowing down exhibited by condensed mat- 
ter systems and reproduced in the theory of dynamic critical 
phenomena Q. 

The failure of (resummed) perturbation theory close to the 
critical point is by no means unexpected. When T ~ Tc the 
effective expansion parameter, i.e. \T/m(T), diverges, and 
it is well known that the - super-daisy resummed - effective 
potential is not even able to reproduce the second order phase 
transition of the real scalar theory Q . 

In refs. Jic|-|l^] it was shown that the key effect which is 
missed by perturbation theory is the dramatic thermal renor- 
malization of the coupling constant, which vanishes in the 
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critical region. This can be understood in the framework of 
the Wilson Renormalization Group. The IR regime of the 
four-dimensional field theory at T = Tc is related to that 
of the three-dimensional theory at T — 0. In particular, the 
three-dimensional running coupling is obtained from the four- 
dimensional one by 

A 3D (A)=A(A)^. 

At the critical point, A3D flows in the IR to the Wilson-Fischer 
fixed point value A^ 7^ 0, so that the four- dimensional cou- 
pling vanishes, 



a(a) -» ±\; D 



(for T ~ T c and A -> 0) . 



The critical exponent governing the vanishing of A is the same 
as that for m(T), so that the ratio X(T)T/m(T) goes to a finite 
value at Tc, and a second order phase transition is correctly 
reproduced. 

We will show that the running of the coupling constant for 
T ~ Tc is crucial also in reproducing the expected critical 
slowing down, i.e. the vanishing of the plasmon damping 
rate. As a first rough ansatz one could just replace the tree 
coupling A with the thermally renormalized one in eq. ([]) 
and readily obtain the expected behaviour. As we will see, 
this gives a wrong answer, due to the logarithmic singularity 
of the on-shell imaginary part when m(T) vanishes. 

In the rest of the letter we will illustrate the computation 
of the plasmon damping rate in the framework of the Wilson 
renormalization group approach for thermal modes (TRG) in- 
troduced in ref. This method allows a computation of 
the damping rate for any value of T, from very high values, 
where perturbation theory works, to the critical region, where 
renormalization group methods like those of j?) are necessary 
to resum infrared divergencies. But the TRG is applicable 
also in the intermediate region, in which none of these meth- 
ods can be applied. 

For a detailed discussion of the formulation of the RG in 
the real-time formalism of high temperature field theory the 
reader is referred to refs. |^,^| • The basic idea is to introduce 
a momentum cut-off in the thermal sector by modifying the 
Bose-Einstein distribution function appearing in the tree-level 
propagator as 



N(ko) JV A (fco) = N(ko)0(\k\ - A) 
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where 9(x) is Heavyside's step function, but smooth func- 
tions can be used as well. Following Polchinski (l4| , the ef- 
fective action and the green functions can then be defined 
non-perturbatively by means of evolution equations in A. For 
A>T the statistical part of the propagators is exponentially 
damped by the Bose-Einstein distribution function, and we 
have the T = renormalized quantum field theory (with all 
the T = quantum fluctuations included). As A is lowered, 
thermal fluctuations with \k\ > A are integrated out. In the 
limit A — > the full thermal field theory in equilibrium is 
recovered. 

We will compute the real and imaginary parts of the (1-1) 
component of the self-energy 



where the real part is the same as that of the self-energy 
appearing in the propagator, Ef 1 (fc;A) = n_R(fc;A) and the 
imaginary part is linked to 11/ in (nj) by IIj = Ed/(1 + 27V) 
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FIG. 1. Schematic representation of the evolution equa- 
tions for the two-and four- point functions. 

The evolution equation for En (A) is represented schemat- 
ically in the upper line of Fig.l. The full dots indicate full, 
A-dependent vertices and propagators, which contain all the 
thermal modes with \k\ > A. The empty dot represents the 
kernel of the evolution equation, which substitutes the full 
propagator in the corresponding leg. Its (1-1) component is 
given by 

#A,n(fc) = - PA (k)s{k )A5(\k\ - A)/V(|fc |) 

where pA(k) is the full spectral function and e(x) = 8(x) — 
0(-x). 

In order to determine the A-dependent four-point function 
we must consider the system containing also the second equa- 
tion in Fig.l. Since the theory has a - possibly spontaneously 
broken - Z2 symmetry, the trilinear couplings are not arbi- 
trary and will be determined from the mass and the quartic 
coupling as we will indicate below. All n-point vertices with 
n>4 will be neglected. 

The initial conditions for the evolution equations are given 
at a scale A = Ao S> T (due to the exponential damping in 
the Bose-Einstein function, A <; 10 T will be enough). Here, 
the effective action of the theory is approximated by 
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where Ma and Aa are the renormalized parameters of the 
T — theory. We are interested to the case in which the 
^-symmetry is broken at T = 0, so we will take [i\ < 0. 

Then, we make the following approximations to the full 
propagator and vertices appearing on the RHS of the evolu- 
tion equations: 

i) the self-energy in the propagator is approximated by a run- 
ning mass, II(fc; A) ~ m\ given by 



m 2 L = ma (if Ma > 0) 



-2/4. (if fj.% < 0) 



En (w ± is, k; A) = En (w, k; A) ± iEu (to, k; A) 



(6) 



it) the four-point function is approximated as 

r i 4) (p0 - -Aa - ir)A (p») , 
Hi) the three-point function is given by 

rf fa) ~ (if 11I > 0) , V-6AaMa (if A < 0) • 

Moreover, all vertices with at least one of the thermal indices 
different from 1 will be neglected. 
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Some comments are in order at this point. In perturbation 
theory, the on-shell imaginary part is given by the two-loop 
setting sun diagram. The evolution equations of Fig.l con- 
tain instead only one-loop integrals, so how can an imaginary 
part emerge? The crucial point here is that the momentum 
dependence of the imaginary part of the four-point function 1 
has to be taken into account (see point ii)). When the latter 
is inserted into the upper equation, IT/ is generated. 

As long as the A-dependent action is in the broken phase, 
trilinear couplings will be present. By neglecting the imagi- 
nary part of the self-energy in the propagators on the RHS 
(point i) above), we loose the contribution to Hi obtained by 
cutting the second diagram in the RHS of the upper equation 
of Fig.l. However, if we restrict ourselves to temperatures 
T > Tc, the trilinear couplings will be different form zero 
only in a limited range of A. They will not contribute in the 
IR, where the dominant contributions to the on-shell imagi- 
nary part emerge (see Fig. 3). The error induced by this ap- 
proximation can be estimated noticing that the contribution 
we are neglecting is of the same nature as the one obtained 
by inserting the last diagram of the second equation in the 
equation for the two-point function. The latter is taken into 
account and its effect on the full imaginary part is of the order 1 " 
of a few percent. 

We have now a system of four evolution equations for 
n^(fc;A) ~ m\, n/(fc;A), Aa, and r)A(Pi), with initial con- 




ditions m A 



-2^i Ao , A Ao , and Hr(fc;A ) = r]A {Pi) = 0, re- 



spectively. Moreover, the subsystem for m A and Aa is closed 
and can be integrated separately. So we proceed as follows. 
Fixing the temperature, we first integrate the subsystem for 
m A and Aa down to A = in order to find the plasmon mass, 
m A= o- Then we fix the external momentum of Hi on this 
mass-shell, k = (m\ =0 ,0) an d integrate the full system from 
A = Ao down to A = 0. 

In Fig. 2 we plot the results for the damping rate 7 and the 
coupling constant at A = 0, as a function of the temperature. 
The dashed line has been obtained by keeping the coupling 
constant fixed (A-independent) to its T = value (A = 10 -2 ), 
and reproduces the divergent behaviour found in perturbation 
theory (eq. (^)). The crucial effect of the running of the cou- 
pling constant is seen in the behaviour of the dot-dashed line, 
which represents the main result of this letter. For tempera- 
tures close enough to Tc, the coupling constant (solid line in 
Fig. 2) is dramatically renormalized and it decreases as 

A A =o(T) ~ t" 

where t = (T — Tc)/Tc and we find v ~ 0.53. The mass also 
vanishes with the same critical index. The decreasing of A 
drives 7 to zero, but with a different scaling law, 

7 A=o(T) ~ i"logt . 

The above expression can be understood noticing that the 
two- loop contribution to Hj, computed at vanishing uj = 
m(T), goes as A 2 log m(T). The RG result replaces A with 
the renormalized coupling, then from (fiT) we have 



t=(T-T c )/T c 



FIG. 2. Temperature dependence of the coupling constant 
(solid line), and of the damping rate with the effect of the 
running of A included (dash-dotted) and excluded (dashed). 
The values for 7 have been multiplied by a factor of 10. 

Taking couplings bigger than the one used in this letter 
(A = 10 -2 ), the deviation from the perturbative regime starts 
to be effective farther from Tc- Defining an effective temper- 
ature as A a =o(T)/Aa < 1/2 for T c < T < T e ff we find that 
t e ff scales roughly as t e ff ~ Aa . 

In Fig. 3 we plot the running of Aa(T) and 1a(T) for two 
different values of the temperature. When T Tc most of 
the running takes place for A <; A so f t , so it is safe to stop 
the running at this scale neglecting the effect of soft loops. 
When T — > Tc this is not possible any longer, since most 
of the running of A and 7 takes place for A £ A so / t . Thus, 
the vanishing of the mass gap forces us to take soft thermal 
momenta into account. 
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FIG. 3. A-running of the coupling constant A and of the 
damping rate 7 for two different values of t = (T — Tc)/Tc- 
The values of 7 (dot-dashed lines) have been multiplied by 
500. 

In conclusion, we have obtained the critical slowing down of 
long wavelength fluctuations in the framework of high temper- 
ature relativistic quantum field theory. The dominant effect 
has been shown to be the thermal renormalization of the cou- 
pling constant, which turns the divergent behaviour of the 
plasmon damping rate into a vanishing one. The notion of 
quasiparticle turns out to be valid much closer to the critical 
regime than the perturbative result would indicate. Indeed, 
for A = 10 -2 , the ratio j/m(T) becomes larger than unity 
for t ^ 10~ 5 in perturbation theory, whereas the RG result is 
still 7/m(T) ~ 0.3 at t ~ 10" 9 . 

In a cosmological setting, the increasing lifetime of the fluc- 
tuations of the order parameter may modify the dynamics of 
second order - or weakly first order - phase transitions. In- 
deed, if the thermalization rate of long wavelength fluctua- 
tions exceeds the expansion rate of the Universe, the phase 
transition will take place out of thermal equilibrium Jl(| . This 
would lead to a scenario for the formation of topological de- 
fects similar to that discussed by Zurek in Jl^| . 
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